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The direct CP asymmetry acp(D° —> KsKs) involves exchange diagrams which are induced 
at tree level in the Standard Model. Since the corresponding topological amplitude Ekk can be 
large, D° —» KsKs is a promising discovery channel for charm CP violation. We estimate the 
penguin annihilation amplitude with a perturbative calculation and extract the exchange amplitude 
Ekk from a global fit to D branching ratios. Our results are further used to predict the size of 
mixing-induced CP violation. We obtain | acp(D° —¥ KsKs )| < 1.1% (95% C.L.). The same bound 
applies to the nonuniversal part of the phase between the D — D mixing and decay amplitudes. If 
future data exceed our predictions, this will point to new physics or an enhancement of the penguin 
annihilation amplitude by QCD dynamics. We briefly discuss the implications of these possibilities 
for other CP asymmetries. 


I. INTRODUCTION 

While direct CP violation (CPV) is well established 
in the down-quark sector BHU, CPV has not yet been 
observed in the decays of up-type quarks. For the discus¬ 
sion of CPV in some singly Cabibbo-suppressed D decay 
it is convenient to decompose the decay amplitude A as 

A = X sd A sd - yA. (1) 

Here X q = V* q V uq and A s d = (A s — A<j)/2 comprise the 
elements V l3 of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix. In the limit Af, = 0 all direct and mixing-induced 
CP asymmetries vanish in the Standard Model (SM). The 
suppression factor Im ~ —6T0 -4 makes the discovery 
of CKM-induced CPV challenging. At the same time 
this parametric suppression renders CP asymmetries in 
charm decays highly sensitive to physics beyond the SM. 

In this paper we study the decay D° —► KsKs . For 
this decay mode A s d vanishes in the limit of exact SU(3 )f 
symmetry HME so that the branching ratio is sup¬ 
pressed. However, Ab does not vanish in this limit and we 
expect \Ab/A s d\ to be large. Therefore CP asymmetries 
in D° K S K S may be enhanced to an observable level, 
even if the Kobayashi-Maskawa phase is the only source 
of CPV in charm decays Hsus]. Moreover, a special 
feature of D° —»• KsKs is the interference of the decays 
cu —> ss and cu —>• dd, both of which involve the tree- 
level exchange of a W boson (exchange topology E , see 
Fig. 0- This interference term gives a contribution to Ah 
owing to Xd + A s = — A;,. That is, contrary to the widely 
studied decays D —► 7r + 7r~, 7r°7r°, K + K~ , no penguin di¬ 
agrams are needed for nonzero direct or mixing-induced 
CP asymmetries. Moreover, the exchange diagram E is 
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enhanced by a large Wilson coefficient. These properties 
make D° —► KsKs an interesting discovery channel for 
CPV in the charm system. 

In this paper we calculate the allowed ranges for the 
direct and mixing-induced CP asymmetries in D° —> 
KsKs, using the results of our global analysis in 
Ref. PE], There are two ingredients which we cannot 
extract from this analysis: the first one is the penguin- 
annihilation amplitude PA (see Fig. [I]), which we esti¬ 
mate with the help of a perturbative calculation. The 
other undetermined quantity is a strong phase 6, whose 
value, however, can be determined from the data once 
both the direct and mixing-induced CP asymmetries are 
measured. The actual size of 6 is not crucial for the 
potential to discover charm CPV in D° —► KsKs'. de¬ 
pending on whether sin <5 is large or small either the 
direct or mixing-induced CP asymmetry will be large. 

Our paper is organized as follows: in Section |n] we de¬ 
rive handy formulae for direct and mixing-induced CP 
asymmetries in terms of A s d and Ab■ In Section |III| we 
relate the CPV observables to topological amplitudes. 
Subsequently we estimate the penguin annihilation con¬ 
tribution, which cannot be extracted from a global fit to 
current data, with perturbative methods in Section |IV| 
In Section [V] we present our phenomenological analysis. 
Finally, we conclude. 


II. PRELIMINARIES 

In this section we collect the formulae for the CP asym¬ 
metries. We write 

A(D° -a K S K S ) = ~^=A(D° ->■ K°K°) , (2) 

accommodating the Bose symmetrization of the two Ks’s 
with the factor of \/\/2. Here we identify K$ = (K° — 
K°)/y/2 and assume that the effects of kaon CPV are 




2 


properly subtracted from CP asymmetries measured in 
D° K s K s , as described in Ref. P33- Adopting the 
convention CP|D°) = — | D°) [TSj the amplitude of 
D° -> K S K S is 

A = -K d A sd + yA. (3) 

The direct CP asymmetry reads 


dir = \A\ 2 - \A\ 2 
CP \A\ 2 + \A\ 2 


Im Ab 

w 


Im yr-|Ad| ■ 

•s^-sd 


(4) 

(5) 


Here and in the following we neglect terms of order X 2 
and higher. Furthermore we use the PDG convention for 
the CKM elements, so that A sd is real and positive up to 
corrections of order Af,. 

For the discussion of mixing-induced CPV we also fol¬ 
low the conventions of Ref. [T8j: with the mass eigen¬ 
states |-Di, 2 ) = p\D°)±q 

(j) governing CPV in the interference between the D — D 
mixing and the D° —» KsKs decay through 




D ) we define the weak phase 


q A 
p A 


vKd 
P x sd 


1 - 


1 - 


Ag A b 

2^*sd Asd 

A b Ab 

2A sd Asd 


q 

A 

p 

A 


( 6 ) 


In this paper we focus on CPV effects which are specific 
to the decay D° —>• KsKs . It is therefore useful to de¬ 
fine a CP phase </> m i x which enters all mixing-induced CP 
asymmetries in a universal way: 


P X sd 

Comparing Eqs. Q and ([7|) one verifies that 0 m i x co¬ 
incides with (j> if one sets A & to zero in A/A. In the 
hunt for new physics (NP) in D — D mixing, which 
may well enhance (f> m i x over the SM expectation 0 m i x = 
0(lm Xb/Xsd), one fits the CPV data of all available 
D° decays to a common phase 0 mix [US US- In the 
case of D° —► KsKs, however, we face the possibility 
that already the SM contributions lead to the situation 
|0| > |0mix|- Comparing Eq. § with Eq. 0 one finds 

^ _ K Ab 

2 Ksd Asd 

1 A& Ab 
2A sd A sd 

= (1 - a^ P ) e^"^), (8) 

where we have used Eq. discarding higher-order 
terms ~ (a^p) 2 as usual. By expanding Eq. (|8j) to first 


(fr mix) 


0 i<t>n 


(7) 


order in A& and 0-0. 

A b 

Xsd 


mix Im , Re 


mlx we arrive at 
Ab 

A sd 


“T“rr" Re-y—|A s d|. (9) 

\A\ A sd 


Eqs. <[5j) and Q form the basis of the analysis presented 
in the following sections. In Eqs. (5j) and ([ 9 ]) |M| is triv¬ 
ially related to the well-measured branching ratio: 



16 ttiti 2 d 


\/"ifl - 4m|- 0 , (10) 


The experimental value is B(D° a- KsKs) = (0.17 ± 
0.04) • 10” 3 [21j. The nontrivial quantities entering the 
predictions of a^p and <j> — cf > m ; x are Ab and the phase of 

A sd - 

The time-dependent CP asymmetry reads 


Acp{t ) = 


T(D°(t) -+• K S K S ) - T(D°(t) -> K S K S ) 
A K S K S ) + T(D°(t) -A K S K S ) 
t 


= a 


^dir 

CP 


-Ar 


Here r is the D n lifetime and 
2 \ 


Ar= 


V-( 

q 

K 

p 


-1 - 4 


(ii) 


y cos 0 — x sin 0. (12) 


Eq. (12) contains the mass difference AM and the width 


difference AT between the mass eigenstates D 1 and Z) 2 
through x = tAM and y = rAr/ 2 . In Eqs. (Ill and 
(12) all quadratic (and higher) terms in tiny quantities 
are neglected. In time-integrated measurements, LHCb 
measures the quantity 


Acp = 0 cp — A v —, 

T 


(13) 


where (t) is the average decay time. CLEO has mea¬ 
sured I2U 


A^f° = -0.23 ±0.19. 


(14) 


Recently LHCb has reported the preliminary result [ESI 
Acp Cb = -0.029 ± 0.052 ± 0.022. (15) 


III. TOPOLOGICAL AMPLITUDES 


The decomposition of A sd and Ab in terms of topolog¬ 
ical amplitudes reads m 


A sd = 


Ei ± E 2 — £3 


— — A s 


V2 

2 E kk +PA 


(16) 


A b = 2E+Ei+E y E3+fM a?) 


(18) 
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Here Ekk = E + E\+E 2 is the combination of exchange 
diagrams appearing in D° —> K + K~. The exchange ( E ) 
and penguin annihilation (PA) diagrams are shown in 
Fig. [T] F/ 12,3 account for first-order SU(3 )f breaking 
in diagrams containing s-quark lines (for their precise 
definition see Table II of Ref. 16]). As in Ref. |2Bj PA q 
denotes the penguin annihilation diagram with quark q 
running in the loop. We use the combinations urn mum 


X S P A s + X d P A d + XbP Ab — 

X sd (PA s - PA d ) + (PA S 


= X sd P Abreak p ^ ^ ' 


PA d - 2 PA b ) 

(19) 

( 20 ) 


We recall that E, Ey 2 ^, PA,... are defined for D° —>• 
K°K° or D° — > K + K~. Since A sd and Ah instead in¬ 
volve KsKs , the factor of — 1/x/2 of Eq. |2J) appears in 
Eqs. (fl 6 ) to (fl 8 |). 


Next we define the strong phase 

(2 E kk + PA\ 
A, d )' 

and the positive quantity 

_ ImAb 12 E K k + PA\ 

T7 71 ' 

With Eq. (18 1 we can write Eq. Q as 


( 21 ) 


( 22 ) 


1 yX b —A sd + (2E KK + PA )/\/2 

Pmix = Im -— Re- 


A sd 

= — Im-^— RcosS . 

Xsd, 


Asd 


In the same way one finds 

Oq t p = —Rsm.5. 


(23) 


(24) 


Eqs. (231 and (24) mean that a^p and </> — </> m i x lie on a 


circle with radius R centered at (—Im^-, 0). The allowed 
points are parametrized by the phase S, which we cannot 
predict. The actual value of 5, however, is of minor im¬ 
portance for the discovery potential of CPV, because S 
only controls how the amount of CPV is shared between 
The crucial parameter is R, which 
determines the maximal values of |a^!p| and \4> — </> m ix|- 

are precisely measured one can 


cr and 4 > - </>n 


Once a^Jp and (j) — <j), 
determine R through 


R = 



- 0mix + Im 


AfW 

Asd / 


2 


(25) 


The experimental value can then be confronted with the 
theoretical estimate presented in the next section. The 
impact of our estimate on a^p and <j> — 0 m i x will be pre¬ 
sented below in Fig. [4j 




FIG. 1. Topological amplitudes: (a) exchange ( E) and (b) 
penguin annihilation (PA). Ref. [29f claims that D° —> 
KsKs is Zweig suppressed, but this statement is only true 
for the PA diagram. 


IV. ESTIMATE OF I PA AND R 


The quantity \Ekk \ can be determined from our global 
fit to branching ratios [XB] . For the calculation of PA 


we exploit the large momentum 
through the penguin loop in Fig. 1(b) 


^ 1.5 GeV flowing 

_ and calculate this 

loop perturbatively as in Ref. m, Such methods are rou¬ 
tinely used in B physics [H(ll-f55] . but their applicability 
to charm physics is not clear. 

We work in a five-flavor theory, so that only current- 
current operators appear in the effective Hamiltonian. 
With Q 2 _= (us)v-a(sc) V -a + (ud) V -A(dc)v-A - 
2 (ub)v-A(bc)v-A and the Wilson coefficient C 2 ~ 1.2 
we may write 

A b = ( ^C 2 (K s Ks\Q 2 \D 0 ), (26) 

because the contribution of the color-flipped operator 
Q 1 is highly suppressed. For our estimate of the ra¬ 
tio PA/Ekk in this section we adopt the SU(3 )f limit 
and identify Ekk with E. In this limit we can combine 
Eqs. (26) and ( fl7| ) into 

°^C 2 (K S K S \Q 2 I D°) = 2E ^[ A ■ (27) 

The penguin diagram can be written as [361 
u s 


PA = G f —C 2 x 
A ir 


6 

x V (4 + 4 - 24) (K S K S \Qi |r>°) 


2=3 


(28) 


with the loop function ?'| 4 = r| 4 (g 2 , m 2 ,4) = r 26 defined 
in Ref. [36!. /i ~ yfq* is the renormalization scale which 
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also enters a s and C 2 in Eq. (28). Qz-q are the usual 


four-quark penguin operators, we will need 

<24,6 = {u a C P )v-A ^2 (q P< f)vTA- (29) 

q—u,d,s,c,b 

PA is color-suppressed w.r.t. E and this suppression is 
encoded in Eq. (28) through a s ~ 1/N C . The contribu¬ 


tions from the matrix elements (< 23 , 5 ) are further sup¬ 
pressed and are neglected in the following. We write 


m 


(Q 6 ) = -2 (My A 


■Ap 


VA, 


with 


Ka = (K S K S | ( q a qp)v(u 0 c a ) A I D°) . (30) 


The other quark flavors in the sum in Eq. (29) contribute 
to D° —► KsKs only through another loop diagram, 
yielding a contribution of higher order in a s . With 


p=r 24 + r | 4 - 2 r 24 , 
we can write PA in a compact form: 


(31) 


PA/E we need My ME. We can place a bound on this 
quantity with Eq. (35), if we assume that jMy^l is not 


much larger than | M^y + My A |; i.e. we do not consider 
the case of large cancellations between M AV and My A 
in E. In view of the fact that E is numerically large [TB] 
this assumption seems justified. Writing 


My A — k(M av 


< a ) 


we vary |«| between 0 and 2. Now Eq. (35) entails 


Ka 


E -G F C 2 (1 + K^ri 4 )' 


and thus 


\2E K k + PA\ = \2E K k I 

< \1Ekk\ x 1-3 


1 °/s_ K 

2n P l + n c ^r^ 


(36) 

(37) 

(38) 

(39) 


PA = ~G F —C 2 pMy A , 

7 r 


(32) 


where we have invoked the SU(3 )_f limit to set My A = 


Mi 


VA . The /x-dependence cancels inp, which furthermore 


does not depend on m c in the considered leading order. 
It is an excellent numerical approximation to expand p to 
first order in m 2 s /q 2 and q 2 /ml (while setting md = 0). 
The expanded expression reads 


10 2 
p = - ITT — 


2m? 


+ 


2 q 2 
15 m 2 


2 q/ 

3 11 m 


(33) 


It is worthwhile to discuss how this result translates into 
an expression in a four-flavor theory, in which the 6 quark 
is integrated out at the scale p b = Q(rrih)' : in this alter¬ 
native approach the piece — 2r 2i of Eq. (28) resides in the 


initial conditions of the penguin coefficients gener¬ 
ated at pb- The four-flavor theory permits the use of the 
renormalization group (RG) to resum the log Ai(p b /\Jq 2 ) 
to all orders in perturbation theory, but this resumma¬ 
tion is inconsistent since In (pb/'\/q 2 ) is smaller than the 
nonlogarithmic terms in — 2r 2i . Without RG summation 
the four-flavor theory reproduces exactly the analytic re¬ 
sult in Eq. (33), which is independent of renormalization 


scale and scheme. 

To estimate My A 


we want to relate it to E using 
Eq. (27). After Fierz-rearranging Q 2 we can express the 
LHS of Eq. (27) in terms of M'L A and 


M\y = (K S K S I (q a qp) A (upc a )v I D°) . (34) 

The exchange topology reads (cf. Eq. ( [27] )) 

E = G f C 2 (K S K S | ( ud) V - A (dc)v-A |-D°) — PA d 
= -G f C 2 (M av + My A ) - G F — C 2 r/; 4 My A . (35) 


To leading order in a s we have therefore E = 
— G F C 2 (M AV + My A ). For the desired estimate of 


Here we have used p = y/q 2 = 1.5 GeV, m s (p ) = 
0.104 GeV, m b (p) = 4.18 GeV, and a s (p) = 0.328. 
(r 24 is evaluated in the NDR scheme.) Inserting finally 
Eq. (39) into Eq. (|22|) gives the upper limit 


R<- 1.3 


ImAf, \2E K k\ 


(40) 


This bound determines the radius of the circle whic h de¬ 
fines the allowed area for i x , a// r P ) viaEq. I.e. 

Eq. (40) determines the maximal size of both |q yip and 
10 — ^mixI (neglecting the small ImAf,/A s d in Eq. (23)). If 
future data violate Eq. ( |40| ), this will signal new physics 
or a dynamical enhancement of PA over the perturba¬ 
tive result in Eq. (32). Sec. [V] discusses how these two 
scenarios can be distinguished with the help of other mea¬ 
surements. 

The relation of rlM 2 
given as 


, rriq, /i 2 ) to G(s, x) in Ref. 


is 


' 24 




2 2\ 1 1 1 ( A 4 

>„«)=3-3 1 og( s? 


-ic 


m q — is 




m* 


(41) 


with an arbitrary mass m 2 . Note that the — ie pre¬ 
scription is essential here; an erroneous omission of this 
small imaginary part results in a numerically large mis¬ 
take. The prefactor of G(x,y) in Eq. (41) disagrees with 
Ref. [14]. We further find that the 6-quark contribution 
—2r 24 is numerically as important as r 24 + rf 4 : 

r 2 4 (9 2 i 0, /Li 2 ) = —0.22 — i 1.05 , (42) 

r 2 4 (« 2 , TO s> P 2 ) = — 0.23 — * 1.05 , (43) 

"~2r 24 (g 2 , m^, ^t 2 ) = —2.02. (44) 


























5 


with 1/JVc 

PA = 0 PA 7^ 0 
\E K k/T U c \ < 1.5 

R < 0.004 0.006 


without 1/JVc 
PA = 0 PA ^ 0 
2.6 

0.009 0.011 


TABLE I. 95% C.L. upper limits (Ax 2 = 3.84), with or with¬ 
out 1/JVc input for C and E. 


V. PHENOMENOLOGY 


The last element needed for the calculation of our 


bound in Eq. (40) is \E K k\- To find \E KK \ we employ 


our global fit to all available branching ratios of D de¬ 
cays to two pseudoscalar mesons m • Note that the main 
constraint on this quantity stems from B(D° —>• K + I\~) 
(see Table III of Ref. PS]). The D decays entering our 
fit involve other topological amplitudes in addition to E 
and PA: in the following we refer to the color-favored 
tree (T), color-suppressed tree (C), annihilation (A) and 
penguin (P) amplitudes. 

We consider two scenarios: in the first scenario the 
SU(3)p’-limit amplitudes C and E are varied completely 
freely. In the second scenario we apply l/iV c counting 
[371439] to the amplitudes, where N c = 3 is the number 
of colors. To leading order in 1/N C one can factorize T 
which results in 

Tfac “ - • ( 45 ) 

Here a\ = 1.06 is the appropriate combination of Wil¬ 
son coefficients, m v and are the mass and the decay 
constant of the pion, respectively, and is the appro¬ 
priate D —>• 7T form factor. (Recall that the SU(3)^-limit 
amplitudes are defined for decays into pions.) In our 
second scenario we assume that |{C + 6A)/T tac \, \(E + 
(5x)/T fac | < 1.3 [26], where 6a parametrizes 1/7V 2 cor¬ 
rections to the factorized annihilation (A) topology [IB] . 

The Ax 2 -profile of \Ekk /T fajC \ returned by our global 
fit is shown in Fig. |2ja) . Fig. [2jb) shows the Ax 2 -profile 
of R for the special case PA~ = 0, in which the whole 
effect comes from the exchange diagram Ekk- The 
corresponding 95% C.L. bounds on \Ekk/T iac \ and R 
inferred from Fig. [2] and Eq. (40) are given in Table [I] 
and illustrated in Fig. [ 3 J Note that we do not treat T fac 
as constant, but also fit the form factor F c P 7r . Likewise 
our fit permits B(D° —> KgKs) to float within the ex¬ 
perimental errors. 

pig. a condenses the main results of this paper into a 
single plot: the radial lines correspond to fixed values of 
the strong phase 6 in Eqs. (23) and (24) in the <j>— 0 m ; x ¬ 


ci^p plane. The red and blue discs show the allowed 
regions for the two considered scenarios. Note that our 
bounds depend on branching ratio measurements only 
and do not involve correlations to other CP asymme¬ 
tries. The black circles correspond to different values 


of \2Ekk + PA\ in Eq. (22). Future data on <j) — j x 
and a^p will allow us to determine 6 and R. The experi¬ 
mental value of R can then be confronted with the upper 
limits in Table [T] to probe the color counting in Ekk and 
our estimate of PA. New physics will mimic a dynamical 
enhancement of PA. In case an anomalously large value 
of R will be found, one can proceed in the following way 
to discriminate between different explanations: 


(i) Several CP asymmetries involve PA, but do not 
grow with |F|. For example, 

a Cp(D° -A K+K~), acp(D° tt+tt"), 
a Cp( D ° 7r ° 7r °), (46) 

all depend on P + PA and are expected to be en¬ 
hanced with PA as well, unless the increase is com¬ 
pensated by — P. But in this case instead 

a^(D+ -> K s K+),af P (D+ -> K s ir + ), 
a^ P (Df ^ K+n°), (47) 


which involve P rather than P + PA , become large. 
Thus a breakdown of color counting in Ekk can be 
distinguished from an enhanced PA. 


(ii) PA can be enhanced by QCD dynamics or by new 
physics. In the first case the CP asymmetries in 
Eqs. (46) and (47) will still obey the sum rules of 
Ref. m New physics will violate these sum rules 
if it couples differently to down and strange quarks. 


We close this section by comparing our result with 
other estimates of a ^ P (D° -4 K S K S ) in the literature. 
Using generic SU(3 )f counting Ref. [13] quotes 


I a d c %(D° -4 K S K S ) \ < ~ 0.6%, (48) 

where e quantifies SU(3 )f breaking. Our result in Table]!] 
agrees with this estimate. However, if the possibility of 
a large, l/iV c -unsuppressed \Ekk\ is realized in nature, 
< 2 p!p can be twice as large. 

Ref. [15] relates a,Qp(D° —> KsKs) to Aa^p = 
a^p(K + K~) — « c ? p ( 7r+7r_ )- With present data this re¬ 
lation reads 


.dir 1 < _ x Aa £r p = 0 4 %. 




(49) 


This estimate assumes that two matrix elements corre¬ 
sponding to different SU(3 )f representations are similar 
in magnitude. We remark that there is no strict correla¬ 
tion between c 1 cp(D 0 —> KgKs ) and Aa^Jp, because the 
two quantities involve different topological amplitudes. 


VI. CONCLUSIONS 

We have studied the direct and mixing-induced CP 
asymmetries in D a —> KsKs in the Standard Model. 
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\Ekk/T [> " \ R with PA = 0 


(a) 


(b) 


FIG. 2. (a) Ay 2 profile of \Ekk\- (b) Ay 2 profile of R (defined in Eq. (22 1 ) for PA = 0. The blue and red curves correspond 
to the scenarios without and with l/N c counting applied to C and E. Note that the red line lies partially on top of the blue 
line. 



PA = 0 PA ^ 0 PA = 0 PA ^ 0 


FIG. 3. Theoretical upper bounds on R. Predictions with 
(without) l/iV c counting are labeled “small E” (“large E”). 
To visualize the contribution from exchange diagrams, we also 
show the result for PA = 0. The case PA 7 ^ 0 is based on the 
estimate in Eq. (401. 


a perturbative calculation and obtain Ekk from a global 
fit to D branching fractions as described in Ref. [16 a . We 
find 


|<4^|<1.1% (95% C.L.), (50) 

10-^mix + Im^-l < 1.1% (95%C.L.). (51) 

A simultaneous measurement of a^'p and (f> — 0 m i x will 
determine |2 E KK + PA\. A violation of the bound 


y a Cp 2 + ~ ^mix + Im ^^) - 

will point to an anomalously enhanced PA. In this case 
other CP asymmetries will also be enhanced. 

Note added in Proof: The authors of Ref. m have 
informed us that they agree with our expression Eq. (ETl) . 
The apparent difference is due to a typo in Eq. (13jof 
Ref. [14] . By comparing our numerical codes we could 
trace our numerical differences back to the “ie problem” 
mentioned at the end of Sec. Hvl 


The allowed region for the corresponding two quantities 
a^p and 0 — </> mlx is a disc whose radius can be calculated 
in terms of the exchange amplitude Ap-p- and the penguin 
annihilation amplitude PA. We estimate PA/Ekk with 
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FIG. 4. Correlation of direct and mixing-induced CPV for D' 
(Ax' 2 = 3.84) upper limits on a'c v P {D° -> KsKs) and f — l: 


KsKs from Eqs. (231 to (25). The one-dimensional 95% C.L. 


are shown in blue. If in addition 1 /N c counting is applied to the 
topological amplitudes C and E, the allowed region shrinks to the red area. The black dashed lines show the radii which are 
obtained when setting |2 Ekk + PA|/2.52 ■ 10 -6 GeV to the annotated values. The chosen reference value is the typical size of a 
factorized tree amplitude in D decays, T fac = 2.52- 10~ 6 GeV. Further, for the black dashed lines B(D° —> KsKs) = 0.17-10 -3 
is used. The circle is centered at (—Im(A6)/A s( j, 0) = (6 • 10~ 4 ,0). 
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